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Abstract 

We have developed a versatile software package for the simulation of di-electron pro- 
duction in pp and dp collisions at SIS energies. Particular attention has been paid to 
incorporate different descriptions of the Dalitz decay A — > Ne^e~ via a common in- 
terface. In addition, suitable parameterizations for the virtual bremsstrahlung process 
NN NNe~^e~ based on one-boson exchange models have been implemented. Such 
simulation tools with high flexibility of the framework are important for the interpre- 
tation of the di-electron data taken with the HADES spectrometer and the design of 
forthcoming experiments. 

1 Introduction 

Experiments with the High Acceptance Di-Electron Spectrometer (HADES) [1] are aimed at 
searching for medium modifications of hadrons at high density and moderate temperatures cre- 
ated in heavy-ion collisions in the 1-2 AGeV impact beam energy range. There exist a multitude 
of predictions, partially conflicting in details, awaiting verification or falsification [2, 3]. Due to 
negligible final-state interactions with nuclear matter, di-electrons or di-muons are considered to 
be useful penetrating probes for this purpose. 

While at higher energies various experimental set-ups for di-lepton measurements have been 
or are operating, e.g., HELlOS-3, CERES/NA35, NA38, NA50, NA60 (the latter three for di- 
muons), and PHENIX, HADES is the only presently active installation in the relativistic regime. 
In addition, it can cover elementary hadron reactions (pp, np, and via tagging a spectator, also 
pn) and hadron-nucleus {pA, nA) collisions. This large range of reactions is related to the 
capabilities of the Schwerlonen-Synchrotron SIS 18 at GSI, Darmstadt. 

When searching for medium modifications of hadrons in the di-electron channel, it is impor- 
tant to have a reliable experimental reference, in particular from elementary hadronic reactions. 
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This became clearly evident in view of the unexplained pair excess measured by the DLS and, 
more recently, HADES experiments [4, 5]. In addition, the knowledge of the elementary pro- 
cess A^A^ NNe^e~ is a prerequisite to understand possible in-medium effects in heavy-ion 
dilepton data [3, 6]. In this context, HADES has performed two di-electron experiments us- 
ing a liquid hydrogen target and proton/deuteron beams [7, 8]: pp at 1.25 GeV and dp at 1.25 
AGeV, i.e., at the same kinetic beam energy per nucleon, which access a broad range of topics. 
For example, the branching ratio and involved electromagnetic transition form factors of the A 
Dalitz decay (A N'j* A^e+e~) are unmeasured. In particular, the di-electron production 
in the A^A^ collision is regarded to be sensitive to the nucleon form factor in the time-like re- 
gion [9, 10]. Moreover, the cross section of non-resonant virtual photon emission (often referred 
to as "bremsstrahlung") differs by up to a factor 4 in the most recent calculations [11, 12, 13]. 

On the other hand, the different contributions of short-lived sources are not easy to separate, 
as, in principle, they have to be treated in a coherent approach, which is done usually in quantum 
mechanical calculations including the interferences, using e.g. One-Boson Exchange (OBE) 
models. Such calculations for the process A^A^ — > NNe^e~ have already been done in [14] and 
more recently in refs. [11, 12, 13]. 

Our goal is to present the methods and their applicability to the HADES pp/pn data in order 
to have a simulation tool at our disposal to be sensitive to additional sources going beyond the 
A Dalitz decay contribution. Due to restricted phase space coverage and efficiency, a flexible 
simulation tool, which is capable to make direct comparison of model predictions with data is 
particularly useful in this respect. Here we describe an extended version of the event generator 
Pluto [15] which is the standard simulation tool for the HADES experiments. 

The calculations mentioned above ensure a coherent treatment of the A^A^ NNe'^e^ reac- 
tion, including graphs involving nucleons, A's or higher resonances and fulfill gauge invariance. 
As the graph involving the A Dalitz decay process is expected to be dominant in the pp reaction 
and still important in the pn reaction, a separate treatment of this contribution is useful. There- 
fore, one has to consider two mechanisms for simulations: either a full calculation including 
properly the interference effects, or the production via resonances (e.g. the A(1232)) and their 
subsequent decay, the so called A Dalitz decay model. 

Following these aims a versatile simulation framework has to be able to: 

1. convert parameterized (or calculated) differential cross sections of the A'^A'^ NNe^e" 
reaction into "exclusive" events and, alternatively, 

2. produce di-electrons via resonances (A^A^ — > AA^ — > NNe^e~) using mass-dependent 
branching ratios and angular distributions. 

For each of these two methods different descriptions should be compared. These new develop- 
ments are useful wherever simulations of this kind (in A^A^ as well as in future -kN experiments) 
have to be done in the context of the interpretation of HADES data. 

The main goal of this report is to describe a standardized method to incorporate calcula- 
tions for NN — > NNe^e' reactions wherever available and subsequently to compare them to 
experimental data using an open-source and adaptable package, and to demonstrate how simu- 
lations of the dp/pp reaction may be performed to provide a valuable tool for interpreting the 
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HADES data [8]. Our paper is organized as follows. In section 2, we outline the numerical 
implementation, the software framework and how the simulation of both types have been done. 
The simulation results are discussed in section 3. Our summary can be found in section 4. In 
appendix A, we explain the models included to describe the A Dalitz decay and the used form 
factors. 



2 Numerical realization 

2.1 Pluto framework 

The simulations which are presented here have been elaborated within the context of the Pluto 
framework [15] originally intended to be used for experiment proposals. The Pluto package is 
entirely based on ROOT [16] and steers the event production with very little overhead by using 
so-called "macros" which are - within the ROOT framework - based on the C++ language. 

After the set-up procedure, the event loop is called N^^ times which creates the momenta of 
all involved particles (and masses in the case of unstable ones) event-by-event. Subsequently, 
each event is usually filtered with the detector acceptance or fed into a full digitization package 
like Geant [17] (not part of our package). The pp (or pn) system enters - in our case - as a 
"seed object" into the decay chain with a given center of momentum (cm.) total energy and 
momentum. Cocktails - which are the incoherent sum of different reaction channels - can be 
generated as well. 

Recently, the Pluto package was re-designed in order to introduce a more modular, object- 
oriented structure, thereby making additions such as new particles, new decays of resonances 
and new algorithms up to modules for entire changes (plug-ins) easily applicable [18]. 



2.2 The A mass shape 

One method of event-based simulations is to set up a reaction of consecutive decays, like pp — > 
pA+ PPlee ~^ ppe^e~ . Hence, this means that the A mass shape (and mass-dependent 
branching ratios) must be known prior to event sampling. 

How does one usually generate A^ev events for the Dalitz decay A — > N^*/! In the first step, 
the A mass shape is sampled, i.e., for each event i, a mass m^^ is assigned. 

Following the usual ansatz (see e.g. [19]) we use the relativistic form of the Breit Wigner 
distribution: 

9 [mA) ^(Ml-mi)2 + mi(rtot(rnA))' 

where ttia denotes the actual energy (resonance mass), and Ma is the static pole mass of the 
resonance. The mass-dependent width is the sum of the partial widths 

N 

r'°\mA) = J2^\mA) (2) 
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with N the number of decay modes. The factor A is chosen such that the integral is statistically 
normalized, / dmA g^irriA) = 1, i-e., eq. (2) leads to the following condition for the mass- 
dependent branching ratio for each decay mode k 

The width for the dominating hadronic decays A Nn is derived from a well-known 
ansatz [19, 20, 21]: 



niA Vg^(^A)/ 



X 



AM a) J 

The dependence on the two decay products with masses rriN and enters via the terms q^{mA) 
and q^{MA), namely the momentum of one out of the two decay products in the rest frame of 
the parent resonance. We follow ref. [19, 21] which uses for the resonance the cutoff parameter- 
ization 

u{mA) = , . , TTTT (5) 

with the parameter /3=300 MeV. 

To simulate the Dalitz decay, the mass-dependent branching ratio has to be taken into account 

as 

Thus, the partial decay width 

p jpA^Ne+e~ 
^A-^Ne+e- N _ / j_ 



r--^^^^-^ (mA) = / dm,,^ (7) 

has to be used in the numerator of eq. (1) which is calculated by integrating eq. (7) with a small 
step size for each tua mass bin, respectively. 

The differential A Dalitz decay width dT/drriee depends on three transition form factors, as 
explained in more detail in the appendix. In particular, we compare the effect of using constant 
form factors, with same values as for a real photon emission ("photon-point" form factors) with a 
two-component model from lachello and Wan [22] which is in line with the vector-meson domi- 
nance (VMD) model. This leads to different mass shapes g^^^<^^<^ (wia) which is exhibited in 
Fig. 1. 

One can see that the mass-dependent branching ratio decreases with mass. This feature is 
an important issue for the di-electron production as large contributions of the higher-mass tail of 
the A resonance and the mass-dependence in the A Dalitz branching ratios affect not only the 
di-electron invariant-mass spectral shape but also the di-electron yield, compared to the hadronic 
channels. Therefore, the mass-dependent branching ratio must be considered even if the effect 
might be suppressed in part by the limited phase space. 

How can an event generator take that behavior into account for the mass sampling? An 
elegant way by using the weighting method is explained in the next subsections. 
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Figure 1: Free spectral shape as a function of ttia of A+(1232) (solid line) compared to 
the distribution functions for dedicated decay states. Short dashed line: A+(1232) e^e^p 
(ref. [24] with "photon point" magnetic form factor), long dashed line: A+(1232) — > 7r° +]? and 
dashed-dotted line: A'*"(1232) 1 + P- In addition, the effect of the two component quark- 
model transition form factor of ref. [22] is indicated by the short dotted line. 



2.3 Weight-based method 

Our first step is to bring the spectrum (represented by a histogram) onto an absolute scale such 
that one can assign each of its bins to a differential cross section. Therefore, total cross section 
models for elementary NN collisions have been implemented. 

For the reactions discussed here, two major contributions (beside the bremsstrahlung) have 
been added, which is the A production [19] (we assume crpp^pA+ = |c^pp^pp,ro) and the close- 
to-threshold T] production [25, 26, 27] (the latter one is needed for the pn case). In order to 
obtain these cross sections independent of the number of events, first a default weight of 1/A^ev 
is applied to the NN seed object. In the decay algorithms, where mass and momentum sampling 
takes place, the weights of all attached distributions are multiplied with the default weight of the 
parent particle to an event weight Wi. In particular, parameterizations for the total cross section 
can be included. This leads - in a 1-step decay as defined above - to the simple relation 



where X is an unstable particle {e.g. 7r°, rf) and Z stands for the remaining particles. Eq. (8) 
is the basic definition to be used in all weighted Pluto simulations: The integral (sum) of the 
resulting histogram weights represents the total (or partial) cross section of an exclusive reaction. 
Therefore, the simulated spectra can be compared directly to the normalized experimental data. 




(8) 
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If we now extend this definition to a two-step process, where a decay X ^ ab follows, the 
weight of particle X has to be folded with the mass-dependent weight of the consecutive decay 
model, which will be described in Sec. 2.6. 

In the simplest case we can consider that production and decay of X are independent, which 
is the case for the rj as well as for the 7r° as they have a comparable small width. In particular, by 
using a sampling model (returning a random event from the known distributions and momenta 
of a, b) this weight is the static branching ratio W-^^"''' — b^^"-^ yielding 



-^ev 



2.4 Flat di-electron generator 

Very often experiments are concerned with regions of the phase space, where a small number of 
events is expected compared to the overall number of a given process. This is clearly the case 
for the electromagnetic Dalitz decays, where the di-lepton yield spans many orders of magnitude 
with a high differential cross section for the low-mass pairs, whereas the high-mass pairs have 
a much lower cross section da/druee- Obviously, a large number of events have to be sampled 
before an acceptable number in the high-mass region has been collected. On the other hand, the 
Monte-Carlo simulations presented here need an adequate statistics for the high-mass region of 
interest. 

The solution used here is that sampling is done using a flat di-lepton distribution first. Then, a 
weight is calculated using the same physics decay model, which was employed for the sampling 
in Sec. 2.2. This means, the decay weight W^^"''^ changes from event to event, depending on 
the values mi = rriee ■ 

By extending eq. (9) we get for the Dalitz decays of the pseudo scalar mesons 



with W^^''"°^^'^~^'^ ^{iTiee) as the differential cross section from [23]. The normalization factor 

^(j7,7r0)^e+e-7 _ ^(j7,7r0)^e+e-7 ^yj^(,,,7r0)^e+e-7 (^l^ 

is used to normalize the spectrum to the selected branching ratio and thereby correct for the fact 
that more events are created in the phase space region with small probability, where the model 
weight is small compared to the region usually containing a large number of events. Here, the 
average decay weight 

TV' =Arev + iVp,e (13) 
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is first calculated for a selected number of events (A^pre = 1000 turned out to be sufficient to 
avoid artefacts) and then adjusted within the running event loop. 

For broad resonances, as discussed in Sec. 2.2, the branching ratio should be an outcome, 
and no precondition of the calculation. The same argument can be brought forward for the 
calculations from refs. [11, 13]. Integrating these distributions is difficult since the underlying 
calculations are only presented for invariant di-electron masses larger then 50-100 MeV/c^ [11, 
13]. Thus, the partial cross section can not be calculated correctly and added to the data base; 
at least it would require some extrapolation. However, this can be avoided, as shown in the 
following. 



2.5 IModels returning dcr / dm 

As an application for a simulation without an explicit re-normalization to a fixed branching ratio 
or total cross section, the above-mentioned weighting has been exploited because the calculations 
from refs. [11, 13] provide the differential cross section da/druee already done on an absolute 
scale. Aiming for a comparison of the A Dalitz decay with the resonant A^^A terms (method 1), 
we define as a model weight W^^{me}) = da/drriee parameterized as described below. It is 
evident that the same method can be used for the full (coherent) differential cross section and 
the quasi-elastic term as well, just be replacing W^^{mel) in the generator by W^'^^^{mel) and 
p^eiaj-^W-j^ respectively. 

In such a case it is more convenient to use the function directly with the flat di-electron 
distribution generator, as in the previous example, but without the intermediate step of production 
and decay. As the event loop generates a row of A^'ev values rrii — rriee the cross section is 
represented naturally by the Monte-Carlo integration method as 



.:^™^+.- f da 



^pp^ppe^e ^ I ^^^^ 

' drUee 

^f^iy^^(m«)-Amee, (14) 



=1 



where Arriee is the kinematic range of the di-electron generator (changing also event-by-event) 
which is provided by the model. It is an advantage that the normalization factor S is not needed. 

The function W^^ required here has been obtained by digitizing the curves provided by 
ref. [11] and using a parameterization of the form 

dm, 



eel 



"ee 



P2 exp(Po"^ee + Pimle 

with 



('^max ^ee) (15) 



mmax = V {mN2 + mNiY + 2mAr27kin - {mm + mN2) (16) 
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as the kinematic limit with A^l as the beam nucleon. The parameters Pi fitted to the curves in 
ref. [11] are polynomials Pi(Tkin) = a° ■ ajTy^in ■ ajT^i^- Th^ function da/dniee of ref. [13] was 
directly supplied by the authors [28] for a fixed kinetic beam energy of Tkin =1.25 GeV, and 
moreover of Tkm =1 GeV and Tkm =1.5 GeV for the pn case. 



2.6 IModels returning dV / dm 

Let us continue with an alternative (method 2), where the A production is followed by the Dalitz 
decay. Pluto treats this process in two steps, leaving out the last and uncritical decay 7*g — > e+e~. 
The weight iy^^^'''=«(mA, ^ee) is now a function of two masses. Taking eqs. (3,7) the mass- 
dependent branching ratio is obtained by 



druee r(t°t)(mA)' 



dm,,^^, r. (17) 



Similar to eq. (14), the di-electron generator represents the Monte-Carlo integration method 
of the model 

jpA-^ATe+e- ( \ 

W ^'=HmA,mee) = ^-1 . (18) 

drUee 

Therefore the mass-dependent branching ratio obtained as 

is already considered, which means that the A mass shape has to be sampled in the first step with 

the pure function g^{m) without using condition (6). This is ensured by the Pluto framework 
automatically. The effect of the mass dependence, as discussed in Sec. 2.2 becomes now visible: 
The fraction of the di-electron events, compared to the hadronic channels, 

^NN^NA^NNe+e- 
K^kin) = ^NN^NA ' (20) 

is in our case with 6(Tkin =1.25 GeV)=4.96 -10"^ significantly larger then the static branching 
ratio, which is b^^^^^ = 4.19 ■ 10~^ (see A.2.1). This means that the mass-dependent branching 
ratio is an important feature which can not be neglected. 



3 Simulation 

3.1 A production 

From the experimental point of view the consideration of known angular distributions is crucial. 
This is in particular true for the emission of the A resonance in the cm. frame which affects 
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Figure 2: (color online) Resonant di-electron production in Np — > Npe^e~ reactions at 
Tkin=1.25 GeV. Left: pp reaction, Right: np reaction. The solid black curve and short-dashed 
(green) curves are the result of the full Pluto simulation using the mass-dependent A width, its 
production cross section and the dV{m/\) / drriee description from [24] with constant transition 
form factors (labeled with "const."), the short-dashed (green) was obtained with the two compo- 
nant form factor [22] (labeled with "2-comp.") in addition. The long-dashed (blue) curve is the 
calculation from [11], whereas the (red) dot-dashed curve is based on ref. [13] and provided by 
ref. [28]. 



the direction of both protons. The impact on the experimentally measured data comes due to 
the fact that the detection of the proton could be a trigger requirement or, explicitly, enters the 
analysis of a (semi-)exclusive channel, e.g. pp Xpe^e^ . Such effects are often integrated out 
in calculations. 

In the case of the A production, we follow the one-pion exchange model of ref. [29], which 
is in excellent agreement with the data and in detail described in ref. [15]. This adds a strong 
forward-backward peaking of the polar angle of the nucleons with respect to the beam axis, an 
effect which is considered in the two-step Dalitz decay simulations. In the OBE simulations vir- 
tual photons are generated with the differential cross sections taken from the coherent calculation 
like a hypothetical "on-mass" shell particle with a given mass niee- Due to missing information 
the momentum and energy is sampled assuming a three body phase space decay of the A^A^ seed 
object sitting at rest in the center-of-mass with m = mNp. Consequently, the virtual photons are 
emitted isotropically in the OBE, but not in the A simulation. 



3.2 Different descriptions of the A Dalitz decay 

After having all needed pieces at our disposal the difference between the OBE calcula- 
tions [11, 13] and the calculations from [24] (either with the photon point form factors or the 
two componant quark model option) can be studied in a quantitative way. Fig. 2 shows the Pluto 
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simulation for the processes pp — > pA+ — > ppe+e (left) and pn — > A'^A''"'^ — > pne^e (right) 

for the two descriptions mentioned above together with the resonant term Np Npe^e~ from 
the OBE calculation [1 1]. The two OBE calculations use the same — A transition form factors, 
but the latter are different than the "photon point" form factors used in the two step Dalitz decay 
model. The effect on the di-electron spectrum is, however, expected to be lower then 15%. 

It is obvious that the disagreement is much larger. The OBE calculation of ref. [1 1] is larger 
by a factor of 2-4 than the production via a "free" A using the mass-dependent branching ra- 
tio. The latter method is a crude approach, taking into account only one graph neglecting in- 
terferences and anti-symmetrization effects. It is, however, frequently used in transport code 
calculations. 

The VMD calculation comes closer to the OBE model, but undershoots it at low masses as 
well. Surprisingly, we do not come to the same conclusion as ref. [13], where the spectrum 

of the two-step Dalitz decay model has a very steep slope. Qualitatively, our model is much 
closer to the OBE calculations of ref. [11, 13] and, taking the VMD form factor into account, 
lies almost on top the curves from ref. [13]. This clearly needs experimental confirmation and 
further theoretical studies. 



3.3 Final state interaction 

In all near-threshold reactions, the final state interactions (FSl) may influence strongly both the 
total cross section as well as the population of the phase space. The first effect is already included 
in the Delta Dalitz two-step simulation as they use measured data. For the OBE calculations the 
factorization [30, 31] 

final ^ yj^fsi ^ yj^AfA,full,ela ^2\) 

has been implemented. This is controlled via a switch which means a factor W^^^ is attached 
exploiting the inverse of the lost function 

/ -I \ 2 / , , . \ 2 



1 y _ f k + ia 
lik)) ~ \k + i(3 



W'^' =1^1 =1 1—^ 1 (22) 

with 

a,f3^-(Jl-'^±l], (23) 



and k as the relative momentum of the two nucleons using the effective radii and scattering 
lengths oo = —7.8098 fm and ro = 2.767 fm for the pp case, and oq = —23.768 fm and 
ro = 2.75 fm for the pn case [31]. We use this description to be compatible with the calculation 
done in [1 1]; other functions could be implemented as well. 



3.4 Nucleon momentum distribution in the deuteron 



For the dp reaction, the nucleon momentum distribution in the deuteron has to be taken into 
account, because the effective neutron momentum may have a big impact at larger di-electron 
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Figure 3: (color online) Differential cross section as a function of the invariant mass rripn in the 
quasi-free pn reaction taking into account the deuteron momentum distribution. Solid curve: 
A+'° production, short dashed (green) curve: pn prni^ir^ (constant cr = 0.1 mb assumed), 
long-dashed (blue) curve: pn pn-q, dot-dashed (red) curve: pn — > dr]. 



masses [12]. This is done by Pluto in a two-step process. In the first step, the off-shell mass of 
the participant 



"^part = + rrtp - 2md ^ m2 + p^^^^ (24) 

is determined by the parameterized wave function poeut from ref. [32]. Along with the second 
reaction particle (in our case the target p at rest) this off-shell particle forms the pn composite 
with a total cm. energy mp„ of the quasi-free reaction. In a second step the bremsstrahlung 
model calculates the energy 

m^„ — ml — ml — 2mpmn 

■^'^ " 2^ ^ ^ 

of the proton in the neutron rest frame. Here, we use the invariant mass mp„ to get the total cross 
section a{mpn) using the on-shell neutron mass m„, which is also the approach in ref. [33]. It is 
important to note that, as the the di-electron cross section for the pn reaction was parameterized 
as a function of the kinetic proton energy, the actual T^in reflects the proton quasi kinetic energy 
even in the case of a deuteron beam and a proton target. 

The consequence of the above-mentioned momentum distribution is that the dp reaction 
results in a "smeared" pn reaction cm. energy and enables thus sub-threshold r] produc- 
tion [25, 26, 27]. The formation of pn composites below any threshold of the final state (as 
e.g. for the pn — > pni] threshold) is rejected while counting the number of rejected events to 
keep the proper normalization. 
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Figure 4: Ratio of (jP"^P"ee^^pp^ppee ^ function of using the OBE calculation [11] with 
the resonant A^A term (dashed curves), the quasi-elastic term (dotted curves) and the full calcula- 
tion (solid curves), for a beam energy of 1.25 GeV/u. The grey (online: red) curves are obtained 
with the pure pn reaction, whereas the black curves include the momentum distribution of the 
deuteron. 



Fig. 3 shows the different channels which are combined later on in the final di-electron cock- 
tails. 

3.5 pn /pp ratio 

In the following, the iso-spin dependence of the bremsstrahlung is discussed. By using the same 
methods for the OBE calculations for the resonant term, the quasi-elastic term and the full cal- 
culation, respectively, and dividing the simulation results for pn and pp the iso-spin dependence 
can be studied. Fig. 4 shows this ratio for the different contributions and the coherent sum, op- 
tionally including the momentum distribution of the deuteron. The simulation indicates that the 
influence of the momentum distribution for masses smaller than 0.4 GeV/c^ is negligible. 

3.6 Cocktail simulation 

In the context of the new HADES data, the simulation has to be done with a full cocktail cal- 
culation. Here, all channels contributing to di-electron production at given energy have to be 
included. The main source of the di-electrons are the 7r° Dalitz decays. The production of vr" 
mesons is done within the A resonance model assuming that all tt'^'s are created via A. The 
production of rj mesons is also included for the dp case. 

Fig. 5 shows the possible contributions for the di-electron production, based on (i) the model 
from refs. [11, 12], and (ii) on the A production and subsequent Dalitz decay. The difference 
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Figure 5: (color online) Cocktail simulation of the differential cross section of da / drriee for the pp 
reaction at 1.25 GeV (left) and the dp reaction at 1.25 GeV/u (right). Long dashed (blue) curve: 
A Dalitz decay using constant transition form factors, short dashed (blue) curve: same with the 
two component quark model from [22, 40], solid black curve: full coherent OBE calculation 
from [11], dot-dashed (red) curve: 7r° Dalitz decay, dotted (green) curve: r] Dalitz decay {dp 
only). In the dp case, the r] contribution haas been added to the short-lived components. 



between these two approaches is clearly visible. Moreover, the Dalitz decay is used with the 
form factor model from lachello and Wan [22]. 

We expect that the HADES data will be sensitive to differences in the various descriptions. 
However, the simulation done with the Pluto framework has to be filtered first with the accep- 
tance of the HADES spectrometer. One limitation is here that the OBE calculations [11, 12, 13] 
provide only the invariant mass spectra, but do not show predictions for angular distributions. In 
the A Dalitz decay model from Pluto, however, these effects are taken into account using existing 
data [15]. We therefore would like to suggest that future theoretical investigations should con- 
sider and present these distributions as this would be extremely helpful for the disentanglement 
as well as understanding of the aforementioned effects, specifically concerning the HADES data. 



4 Summary 

In summary we have presented the details of an extension of the previous Pluto framework able 
to incorporate important descriptions for the production of low-mass di-electrons in elementary 
collisions such as pp and quasi-free pn. Several models based on a free A Dalitz decay and full 
quantum mechanical calculations have been included so far and are ready for the comparison 
with the upcoming HADES data [8]. As the intermediate excitation of A is important for the 
di-electron production, we have described in detail the corresponding models. 
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A A Dalitz decay models 

A.l Differential decay width of the A Dalitz decay 

In this report, for computing the A(1232) Dalitz decay we use the prescription of Ref. [34] 

-lA— >Ne+e 



with a = 1/137 as the fine structure constant. The decay process A N-^f* consists in three in- 
dependent amplitudes, which can be calculated unambiguously from the electromagnetic vertex. 
However, as stressed in [24], inconsistent formula for the differential decay width of this process 
can be found in the literature. We tackled the calculation by ourselves, using the magnetic dipole, 
electric quadrupole and Coulomb quadrupole covariants [35], and could confirm the expression 
of [24], as repeated below: 

a (niA + ruNY i 7 

— ^2 — \ y+y ' 

lb m\mj^ V ^ ~ 
y± = (^A ± rnNf - ml^, (27) 

where the index N refers to the produced nucleon, e is the electron charge, and G^^^* (m-y* ) 
depends on the N- A electromagnetic transition form factors as 

2 

Tfl 

' ' 'G^KOl (28) 



"2< 1^^^"^^*' 

where Gm("^7*)' G£;(m^*), GciiTiy* ) are the magnetic, electric and Coulomb — A transition 
form factors, respectively, which will be discussed in the next section of the appendix. Note 
that eq. (27) implies a normalization of the form factors as in [35], since isospin factors are 
included in the numerical factors. We could also check the validity of the expressions derived 
in [36], where the amplitudes are calculated with a different, but equivalent set of covariants, 
with corresponding form factors. We use eqn. (27) throughout. 
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A.2 Electromagnetic AT — A transition form factors 

The electromagnetic — A transition form factors are analytical functions of the squared four- 
momentum transfers at the — A vertex. Pion electroproduction and photo production ex- 
periments allow to determine these formfactors in the space-like region (g^ <0) [37]. In the A 
Dalitz decay process, due to the positive four-momentum transfer squared (g^ = m^, > 'imlJ, 
the time-like region is probed, where only the limit at = is known experimentally. An addi- 
tional difficulty hails from the fact that the form factors, which are real in the space-like region, 
get an imaginary part in the time-like region. 

Therefore, two options can be chosen to compensate the lack of experimental information on 
these observables and have been implemented for Pluto simulations. 

A.2.1 Constant N — A transition form factors 

In this option, which is based on the smallness of the squared four-momentum transfer in 
the A Dalitz decay process, the form-factors are given the values Gm=3.0, Ge=0, Gc=0. This 
choice is consistent with the precise measurements of the electric and magnetic form factors 
in pion photoproduction experiments [37] and with the very small contribution of the Coulomb 
term in eq. (28) and provides in addition the correct radiative decay width r'^^^'"'* — 0.66 MeV. 
The resulting branching ratio at the pole mass of h^-^^^^ = 4.19 • 10~^ is remarkably consistent 
to the photon decay branching ratio times the fine structure constant a which would result in 

^A-.JVee ^ 4 01 . 10-5. 

A.2.2 Two-component quark model 

The alternative is to use a model for the A?^ — A transition form factors. This model should 
in principle satisfy the analyticity properties in the complex plane, as well as the asymptotic 
behavior predicted by QCD sum rules, while reproducing the existing data measured in the space- 
like region. At first, the photon-point value provides the normalization of the whole function. As 
an example of such models, the two-component quark model, described in the following, has 
been implemented in the Pluto event generator. 

The picture behind this transition form-factor model is the overlay of an intrinsic q^ structure 
and a meson cloud which couples to the virtual photon via vector mesons [38]. The model allows 
an analytical continuation of the function from space-like to time-like region [39], inducing a 
phase, and could be successfully applied to the description of elastic space-like and time-like 
nucleon form-factor measurements. The formalism has been recently extended to calculate also 
baryonic transition form factors in a unified way [22], these new developments have been tested 
on the space-like A^ — A transition form factors measurements. We use here the simplest version 
of the model, which assumes isospin symmetry, and therefore considers only the magnetic A^ — A 
transition form factor. The intrinsic 3 -quark structure is described as: 
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H) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



m^, [GeV/c'] 

Figure 6: Form factor model taken from [22, 40]. 



and the overall expression for the time-like N — A transition form factor is: 

X (P' + PF,{q^)), (30) 

where /Xp = 2.793 is the proton magnetic moment and P and are the constants for the coupling 
to the quark core and to the meson-cloud respectively. In the case of the A^— A transition, only the 
p meson contributes to the latter contribution, due to isospin conservation and the corresponding 
-dependence is given hyFp{q'^), as: 



m'l — q'^ + 4m7r(l — x)Tp{a{x) — i'y{x)) 



(31) 



where we have introduced x — q'^/Aml. and [40] 



a {x) 



2 x-1 



TT V X 



In [y/x - 1 + y/x) 



> if X > 1 



(32) 
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Figure 7: (color online) The distribution 
masses. Solid curves: pole mass m-A = 1.232 GeV, long dashed curves: ttia = 1.5 GeV and 
short dashed curves: rriA = 1-8 GeV. The description from [24] is applied. The upper set of 
black curves has been calculated with the two-componant quark model from ref. [22] whereas 
the lower set of grey (online: green) have been obtained with the constant transition amplitudes 
as described in the text. 



a (x) 



X 



X 



1 cot 

TC 



-1 



X 



1 — X 



> if a; < 1. 



(33) 



7 (x) = J 

In eq. 31, the values mp=765 MeV and rp=112 MeV are used different from the physical 
values due to the form of parameterization. The value of the parameters a=0.29 GeV^^, ^^=53°, 
/?=1.2147 and /5'=0.004 results from fits of the model predictions to the available experimental 
information as discussed above. The resulting distribution of form-factor values (see fig. 6) 
shows a broad peak centered around = rriee ~ 0.6m^. Due to the small value of the (3' 
parameter, the contribution of the coupling to the quark core in this model is negligible up to 
5 (GeV/c)^, the dominant feature of the model in the kinematic range probed by the Dalitz decay 
process is therefore the vector dominance. 

Fig. 7 finally shows the distribution dT^'^^^^'^ /druee for the two form factor models. The 
results clearly exhibit a rising decay width for larger A masses. In a proton-proton collision 
at 1.25 GeV incident energy, the mass of the produced baryonic resonance is limited to 1.48 
GeV/c^, the latter effect will nevertheless affects the shape of the di-electron mass spectrum 
shown in Sec. 3.2. The A Dalitz decay branching ratio, defined at the pole mass is however 
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mainly determined by the values of the form factor at very low q^. The branching ratio obtained 
using the two-component quark model form factor is about 10% larger then with the constant 
"photon-point" value, as discussed in A.2.1. This derives from the fact that the parameters of the 
model are fitted to a set of data over a large range, which results in a slightly too high value 
for the magnetic form factor at — 0. 

One should note that, in the semi-classical description, this would correspond to the pro- 
duction of an on-shell p meson, which should therefore not be added then as an independent 
contribution when enabling the two-component quark model. 
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